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ABSTRACT 



In a wedge shaped ocean, the method ot images is used to 
develop an analytical approximation ot the acoustic pressure 
field. Contemporary work develops acoustic doublets from a 
combination of the source and surface reflection image using 
simple dipole theory. The method of images is then used to 
sum the dipole images. This thesis matches dipole pairs to 
achieve a quadruplet expansion. A computer program using the 
derived quadruplet equation is then created to verity the 
results by comparing them with the "URTEXT" program. 
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I. INTRODUCTION 

This thesis is a continuation in the examination or the 
acoustic pressure field in a wedge shaped ocean. As the 
focus of naval operations shifts to the littoral regions, 
more attention is being paid to the problem ot A.b.W. in 
acoustically shallow water, i.e. regions where the sound 
paths have multiple interactions with the surface and bottom. 

To date, the simplest propagation model has been the method 
of images . 

In its present form, this model assumes isospeed water, 
a pressure release upper surface, and a lossy, penetrahie 
bottom. For parallel surfaces, the source and its images 
form a vertical array with the actual source at the center. 
Furthermore, it tne source is near the surface it can he 
combined with its negative (180° out of phase) surface 
reflection image to form a dipole source. The column ot 
images can then be considered a column or doublet images. 

The equations used in this thesis are the far field 
approximations of the acoustic doublet. An unbalanced 
doublet is defined as having two sources ot nearly equal 
amplitude but opposite phase separated by a distance d. It 
the amplitudes differ greatly, or it the phase difference is 
not 180°; then the doublet approximation is not valid. 

When the method of images is used in a wedge shaped 
ocean, the vertical column of doublet images becomes a circle 
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centered on the apex. In previous work the upper and lower 
image doublets were summed individually with great success. 
This thesis attempts to group upper and lower douDlets into 
double doublets or quadruplets. 

This method can only be used tor extremely small wedge 
angles. With larger angles the sound paths from the upper 
and lower doublet images will be to dissimilar, and tne 
quadruplet approximation will be invalid. Since typical 
slopes for most of the worlds continental snelves are about 
3°, a good analogy would be that ot using a ship mounted 
active sonar searching for a diesel submarine in the shelt 
region . 
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II. WAVE NUMBERS AND SCALING 



Sound propagation in a shallow channel with a pressure 
release surface parallel to a rigid bottom has a propagation 
wave number for each mode defined by the following equation 






71 1 

H fid) 



Where H is the depth of the channel . For a fast 
is the critical angle found from 



U) 

bottom, e 




( 2 ) 



An analogous reference angle 6^ for a slow bottom is 



tan0^= 



( j) 



If the depth H is replaced by the scaling distance. 



R 



R 

X 



( 4 ) 



and the bottom is tilted at angle ; tnen eq . 1 for the 
cutoff value ot the lowest mode witn a fast bottom becomes 



5 

2sin0tanp 

and for a slow bottom 



5 

2tan0tanp 
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III. ACOUSTIC DOUBLET FORMATION 



The sound propagation for spherical spreading of a point 
source is 

^ — ® (7) 

When two point sources of opposite phase are combined to form 
a doublet, equation 7 can be expanded. 

Where the "+" and subscripts refer to the upper and lower 

sources respectively. These values can be expanded 

2\=r+Ar A^^A+AA 

{ y f 

r_=r-Ar A_=Aq~LA 

Where the delta values indicate the incremental ditterence 
tor each source compared to a theoretical point source 
located exactly between them. The pressure equation can then 
be rearranged to yield 

e ( if ) 

^ i-A£ i+A£ 

r r 

An acoustic source located near the surface has a 
surface reflection image of equal amplitude and opposite 
phase. In the far field the image pair is treated as a 
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single aouniet source. 



For a oaianced source, the pressure 



equation can be approximated by 

[g-jXr-Ar) _^-jA:{r+Ar) j 



Ul) 



As seen in fig. 2, Ar is 



Ar= — sin0 
2 



(12) 



or 

Ar^riSinysino 

(13) 

combining eqs . 11 and 13, yields the pressure 

P=2j-^sin {kr.^si nysin0) 1 1 4 ; 

To develop an equation tor the unbalanced doublet, the 
amplitude relationships expressed in eq . 9 must be rewritten 
as 






aa- 



A,-A^ 






So that for an unbalanced doublet 



P=2-— [sin(icr,sin0) -jA-^cos (Jcr,sin0) ] 
r ^ 



U6; 



The distance from the center ot the doublet to the 
receiver can he calculated using the law of cosines. 

6 



T -2 Tj^r 2COS (2j3P±6) 



U7; 



where the ranges r^ and r 2 are measured from the apex to the 
source, and from the apex to the receiver respectively. The 
angle 6 is measured from the surface down to the receiver. 
The range from each successive doublet in a wedge of angle (3 

IS 



Where the angle in the cosine term refers to the angle or 
each reflected image; those with a are the upper images, 

and those with a are the lower images. This equation can 

be rearranged to yiela 



The first order Taylor’s series expansion can then be used to 
simplify eq . i9 



( 2 np* 6 ) 



(lb) 




(ly) 



(ri-r2){l + — [i-cos (2np±8) ] ) 
^2 



1-^2 



(20) 
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IV. QUADRUPLET EXPANSION 



The upper and lower image doublets can be combined to rorm 
a quadruplet. The complimentary doublets are 180° out or phase, 
and the pressure of the new quadruplet is 

Pn--^^(sin[JcriYsin(2np+5) ] 

-sin [icriYsin {2nP~b) ] 

In this case, both of the interior exponents involving Ar are 
negative. The Ar terms are derived trom a comparison to tne 
range trom the primary douoiet, tnereby ensuring that all pnase 
angles are calculated with respect to a common rererence point. 
To find the Ar values used in eq . 21, an approximation or eq . 
gives 

Arr-rJ"- (r 2 ~r, ) - — [1-cos (2 i3Pjl5) ] (2z) 

The cosine term can be expanded and approximatea ror small 
va 1 ues ot 6 . 

COS (2iip±6) -cos (2iip) ±6sin (2rjp) Wd ; 

Equation 21 can now be expanded into the unbalanced quadruplet 
equation . 



(sin [/cr^ysin {2T3p+5) ] c 
-sin [icr^ysin (2np-5) ] e 



U4; 
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It IS convenient to define 



4 )=icriH ( 1 -cos (2/2P) ) 
d=icr^ 6 sin (2/2P) 



So that eq . 24 becomes 

= AZ e 

[sin [icr^ysin { 2 np + 6 ) ] -7-^ cos [icr,ysin ( 2 n|i+d ) J ] e ( 2 t > 

-B^ [sin [icr^Ysin (2np-6) ] -j-^cos [icr^ysin (2nP -5 ) ] ] 

B 

The values tor Aa and Ab must be determined rrom tne reriection 
coefficients . 
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V. 



REFLECTION COEFFICIENTS FOR A SLOW BOTTOM 



The reflection coefficients used in this thesis are 
derived from the Rayleigh reflection coefficient 

i>c-sin0t/sin0_f 

jbc+sin0j-/sin0j 



Where 




(28; 



The angles are the grazing angle 6^ and angle or transmission 
e^. Therefore, 0^ is defined as 










(28; 



Using trigonometric rei at lonsmps , eg. zi can oe 
rewritten . 

bc-Jl-c^cos^Q. / s in0 ,• 

n — ^ 

bc+^h -r^cos^0j/sin0^ 

In addition, the term under the radical can be rewritten. 

-c^cos^0_f-Cy/sin^0^-sin^0^ (^1) 

Combining eqs . 30 and 31, and rearranging yieios 



b-^1- [sin0^/sin0j'^ ^ b-yji-i/x^ 

[sin0p/sinOj] ^ jb+y^l-l/x^ 
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Where x is 






V—fT 

Sin0^ 

For a slow hottom with a very small angle, eg. 32 can oe 
approximated by 

2bc a 
-p=Z^i 

/?=e 

The argument ot the exponential is 

2bc . 2Pz/Pi 






vl-c-* yvCj/c'iK-l 



(3b) 



So that this approximation ot the Rayleigh reflection 
coefficient can be written in a more concise form 

i?=e 

Since each ray intersects the bottom and surtace multiple 
times, a product of reflection coefficients aiong each 
individual path is required. This cumulative coefficient 
will take the rorm 



n (2/3-/n) p] 

. . . 






and the product yields 






( 3 «) 



For very small angles, this can be approximated by 
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(3yj 



The angles of incidence are needed to calculate tne proauct . 
The angles for the bottom images are 

0o-2(ji i)P (40) 

Using the law of sines, the angle at the receiver is defined 

e= (ii/Tq) sin(2/iP + 6) ( 4i; 

The angle of incidence of tne ray from tne center or tne 
doublet image to tbe receiver intersecting tbe apparent 
bottom is 

0=2 (n-l) p-^6 +e ( 42 ; 

However, tne upper and lower images in eacn doublet nave 
sligntly different angles ot incidence witn the apparent 
bottom 

e*=(2ii-l) p-^S-^c, (4^; 

The upper and lower angles at the receiver ror each aouOiet 
are 



t^=e± {i\/i'^)ycos (2iip+6) (44) 

The individual angles of incidence are 

0i=0o+6 (1+ (Xi/Xj,) cos (2i2P + 6) ) ± (Xi/Xo) YCOS (2jjP) ( 45} 

If the following variables are used in the coefficients as 
approximations ot the sine or the angles ot incidence 
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(4b; 



+ ~ ccs2n^ 
b--~ coo2n^ 



then eq. 38 becomes 






(47) 



For each doublet image there is an upper and lower image. 
The upper doublet is detined by with tbe upper image A. 
and the lower image A^ . Likewise, the lower doublet has 
the upper image B. and the lower image B. . 









A. and Aa detined by eq . 15 can be combined with eq. 4« 



A^=Sl„cosh inb„ya) e 



na^ia 

AA=^„siiih {nb„y a) 

When combined, they yield 

-^=tanh {nb^fa) 



B„=8l.cosh (/ub^ya ) 
AB=«„siiili inb„ya) 



-^^=tanh {nb^ya) 

■Or 



(49> 



(bu) 



The result is the full quadruplet equation 






P=Py*Y^-^Kcoshh - 

n«l ^ 

[g-a-j^dgin (jcr^Ysin (2nP+6) ) -jtanhb cos (kTiYsin (2nP+6) ) 
-e«*i^<^gin(icriYsin(2/2p-5) ) -Jtanhb cos {ki^^y sin (2/ip-6) ) ] 
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RESULTS 



VI . 

Three runs at different ranges were made to compare tne 
output ot tne quadruplet model with tnat or tne “URTEXT" 
program. For each run, the source angle is 1°, and the wedge 
angle is 3®. The ratio ot the sound speed in seawater to 
that of the bottom is 1,01695, and the ratio ot the density 
of seawater to that of the bottom is 0.7. This is equivalent 
to a slow bottom consisting ot a combination ot clay and 
ooze; conditions very common in continental sheit regions. 

The runs are tabulated at 0.30° increments ot receiver angle 
from the surface downwards. The results are snown in tables 
1 through 3. At a receiver angle ot zero all three runs 
result in an amplitude of zero, as they sbould. 

At the shortest range. Run 3, the results ot tbe 
quadruplet model differ signi 1 1 cant 1 y trom those or "UKTEXT . 
The percentage dirrerence in amplitude is consistently aoove 
10%. With the exception of the 3° receiver angle, the 
average percentage difference in amplitude tor the other two 
runs is well below 10%. The tar tield doublet approximation 
improves at longer ranges, and the approximations tor Ar 
values become more accurate. Also, the difference in phase 
angle is much greater at the shortest range; half the phase 
angles in Run 3 differ by greater than x/4, but in the other 
two runs, only the phase angle at a 3° receiver angle ditters 
by more than that amount. 



Each ot these data sets is accompanied hy sample plots 
of the quadruplet pressure amplitude and tne reriection 
coefficient for each quadruplet. In each case, tne plots 
were taken with a receiver angle of 0.90°. The magnitude ot 
the complex quadruplet pressure drops oft very quickly with 
the quadruplet number. It is negligible by the bth 
quadruplet. This is a result of the effect ot the rapia 
(exponential) decay of the Rayleigh reflection coetticient 
with increased angle of incidence ot the higher numoer ot 
quadrupl et . 

In practical terms, this means that the quadruplet 
expansion becomes less accurate with larger receiver and 
wedge angles. As these angles increase, the ditrerence in 
the angles ot incidence oetween the upper doublet souna path 
and the lower doublet souna path wiii increase 
proportionately. When these ditterences are raisea to the 
exponential power in the reiiection coet r i ci ent s , tne 
differences in amplituae oetween upper ana lower douoiets at 
the receiver will be quite large. However, the tundamentai 
assumption in the quadruplet expansion is that the upper and 
lower doublet images are equivalent in amplituae, but 
opposite in phase. Therefore, the quadruplet expansion is 
only valid tor very small angles and long ranges. 
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TABLE 1 



BETA=3,GAMMA = 1,R1 = 1,R2 = 100, DENSITY RAT 1 0= . 9 , SPEED RAT X 0= 1 . 0 i t> y 5 
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APPENDIX A 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%QUADRUPLET EXPANSION 

%pquad .m 

%MATLAB 

%Quadruplet expansion ot tne acoustic 
%pressure field in a wedge shaped ocean. 

%Michael Joyce 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

B=input ( ‘ wedge angle in degrees=’;; 

G=input (' source angle from surface in degrees= ‘ ) 
D=input (' receiver angle from surface in degrees= 
rho= input (' water to bottom density ratio='); 
cc=input (' water to bottom sound speed ratio='); 
rl=input (' range of source from apex= ' ) ; 
r2=input (' range of receiver from apex= ' ; ; 

Nl=f ix( 90/B) ; 

B=B*pi/180 ; 

G=G*pi/180; 

D=D*pi/ IbO ; 

%%%%%%%%%%%%%%%%%%%%%%%^%%%%%%%%%% 



%CALCULATE SCALING FACTOR 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

tb=tan(B) ; 
if cc<l, 

tl=acos ( cc ) ; 
t 2=sin( tl ) ; 

else 

tl=acos ( 1/ cc ) ; 
t2=tan(tl) ; 

end 

kl=2*t2*tb; 
k=pi/ kl ; 

%CALCULATE CONSTANTS AND PRIMARY DOUfaLET 

al = 2’*'( 1/rho)/ (sqrt(cc^2-l)); 
r=sqrt(rl 2+r2'2-2*rl^r2*cos(D;); 
r3=abs ( r2-rl ) ; 
r4 = rl’^r2/r3 2; 
mu= r2/abs ( r2- rl ) ; 
q = k*rl*G’^D; 
pl=(2/r)*sin(q) ; 
f = 0; 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 



%QUADRUPLET SUMMATION 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

for n=l : 1 : N1 , 
s ( n ) = n ; 
th(n)=2*n*B; 
d(n)=k*rl*D*sin(th(n)) ; 
phi (n)=K*rl*mu* (i-cos(tn(n)j); 

%%%%%%%%%%%%%^%%%%%%%%%^% 

^KEELECTiUN COEFFICIENTS 
%%%%%%%%%%%%%%%%%%%%%^%^^ 

ru(n) = sqrt ( rl 2 + r2 2“2’'rl’'r2*cos( th( n)+D) ) 

rl(n) = sqrt(rl 2 + r2"2-2*rl*r2*cos(th(n)’-D) j 

r5 (n)=rl/ ru(n) ; 

r6(n) = rl/ rl (n; ; 

ai(n)=l+ro(n;*cos(th(nj ) ; 

bl(n)=r6(n)*cos(th(n) ); 

rr (n) = -B’^al (n 2 ) ; 

R(n)=-exp(rr(n) ) ; 
a(n) = n*al*D’^al(n) ; 
b(n)=n*al*G*bl(n) ; 
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%PRESSURE 



ul(n)=cosh(b(n) ) ; 

u2(n) = exp(“D’^phi(n) ) ; 

u3(n) = (-2/r)*R(n)’^ul(n)*u2(n; ; 



%%%%%%%%%%%%%%%%%%%%%%%%% 

%UPPER IMAGE 

vlln)=sin(th(n)+D) ; 
v2(n) = K’^rl’*^G’^vlin) ; 
v3(n)=sin(v2(n)); 
v4(n}=cos(v2(n} ) ; 
v5(n) = -a(n)“ D’^mu*d( n) ; 
vb(n)=exp( v5tn) ) ; 
v7 (n)=tanh( b(n) ) ; 

vbin) = v6(n)’*^(v3(n)“D’^v7(n)’*^v4(n) ) ; 

%%%%%%%%%%%%%%%%%%%%%%%%% 

%LOWER IMAGE 

%%%%%%%%%%%%%%%%%%%%%%%%% 

wlln)=sinithln)-D); 
w2(n) = k’*'rl*G’^wlUi) ; 
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w3(n)=sin(w2^n) ) ; 
w4(n)=cos(w2(n; ) ; 
w5 ( n) = a( n j + j*mu*d( n ) ; 
w6 (n; =exp( w5 (n ) ) ; 

w7(n)=w6(n)*(w3(n)-D*v7(n)’^w4(n)); 

%%%%%%%%%%%%%%%%%%%%%%%%% 

%SUMMATION 

p(n)=u3(n)’*^lv8(n)-w7(n; ) ; 
t=t+p(n) ; 
pr (n)=real (p(n) ) ; 
pim(n} = imag(p(n ; ) ; 
pz(n)=sqrt(pr(n) 2+pim(nj 2); 
rrru, n ) = - 1 n ; ; 

end 

piot(s,pzj ;gria; 

xl abel ( ' Quadrupl et Number'); 

y 1 abel ( 'Magnatude or Quadruplet Pressure') 

title( 'Quadrupl et Pressure'); 

pause ; 

plot (s , rm) ; grid; 

title( 'Rayleigh Reflection Coet ricient ) ; 
XI abel (' Quadrupl et Number’); 
yiaoeii Coetticient' ); 



31 



pause ; 



ps=pl+ t ; 

ph=imag( ps ) / real ( ps ; ; 
phl=atan( ph ) ; 
phase=180*phl/pi 
f = f ; 

pl=pl; 
ps=ps ; 

P=sqr t ( rea 1 I, ps ; 2 + imagips; Z) 
end ; 
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